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Confinement: Results and Perspectives 

Adriano Di Giacomo®^MCSD]Universita di Pisa and INFN, 
Dipartimento di Fisica, Via Buonarroti 2, 56125 Pisa, Italy*, 


An updated review is presented of lattice investigations of colour confinement. 


1. Introduction 

The existing experimental evidence for confinement is based on the search of fractionally 
charged particles (quarks). Such particles have not been observed, and upper limits have 
been established which are typically 10“^® times smaller than what one would expect in 
the absence of confinement |l| . The natural interpretation of that is that quarks do not 
exist as free particles,and that this absolute property is based on a symmetry[^. 

No direct evidence exists instead for confinement of other colored particles with zero or 
integer charge, e.g. gluons. The question is then: what symmetry of QCD prevents the 
observation of free quarks and possibly of other colored particles? 

A deconfinement transition is observed in numerical simulations of quenched QCD on 
the lattice^]. 

The theory at finite temperature T is described by Euclidean Feynman integral ex¬ 
tending from time 0 to time 1/T, with periodic boundary conditions in time for bosons, 
antiperiodic for fermions. On a lattice AV x iV| this means T = 1/aAV. is the time 
extension. Ns the space extension of the lattice, and Ns 3> Nt- 

In physical units a depends on the bare coupling constant g, at sufficiently small g, as 

1 2N 

a = — exp(-6o/3) P = — (1) 

Al g^ 

Eq.(2) stems from renormalization group and —6o is the coefficient of the lowest term in 
the perturbative expansion of the beta function. Because of asymptotic freedom 6o > 0 
and T = ^ exp(6o/9)- Low temperature corresponds to strong coupling regime (disorder), 
high temperature to weak coupling (order). The symmetry responsible for confinement is 
a symmetry of the disordered phase. 

The way to observe confinement is to look at the spatial correlator of the Polyakov lines 


G(T- y) = {L{x)L{y)) 


( 2 ) 


The Polyakov line L{x) = V exp 
lattice along the time axis. The s 


i Ao(x, t) dr is the parallel transport through the 
:atic potential R(T) between a qq pair at distance \z\ is 
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given by 

one has by cluster property 

G(z) ^ Cexp{-a\z\) + \{L)\^ (4) 

|i]—>oo 

Lattice data show that a temperature Tc exists, such that at T < Tc \ {L) | = 0 so that 
V{z) ~ a\^ (T = string tension (5) 

a linear potential which implies conhnement. 

At T >Tc I (L) I 7 ^ 0 and V(z) ~| 2 i^oo 0. 

The resulting indication is that | (L) | is an order parameter and a a disorder parameter 
for the deconhning transition. 

For SU{2) the transition is 2nd order and belongs to the universality class of the 3d 
Ising model. For SU{3) it is weak hrst order. 

The underlying symmetry is Zn, the centre of the gauge group. In full QCD,however, 
Zat is explicitely broken by the very presence of quarks, and also a does not exist due to 
the so called string breaking, i.e. to formation of gg pairs when trying to pull gg apart. 
A chiral phase transition exists in this case, from the Goldstone phase to the Wigner 
phase and the order parameter is {'ip'ip), the chiral condensate. The connection of chiral 
transition to deconhnement is not clear. 

This situation is unsatisfactory, especially in view of the philosophy of Nq —^ oo. 
According to it the basic structure of the theory is determined at Nc = oo, corrections 
in l/Nc being small and convergent. Quark loops are non leading in this expansion. 
Therefore a good order (disorder) parameter should not depend on the presence of quarks. 

2. Duality. 

To understand the symmetry of a disordered phase the concept of duality is fundamen¬ 
tally. Duality applies to d—dimensional systems, having non local excitations with non 
trivial topology in d — 1 dimensions. 

Such systems admit two complementary descriptions. 

A direct description interms of local helds <F, in which topological excitations p are non 
local. The direct description is convenient in the weak coupling regime (ordered phase). 

A dual description, in which the topological excitations /i are local operators and the 
helds <F non local excitations. The dual description is convenient in the strong coupling 
regime g 3> 1 (disordered phase). The effective coupling of the dual go is related to g as 
go ~ 1/g -C 1. Duality maps the strong coupling regime of the direct description into 
the weak coupling of the dual. 

Duality applies to a number of statistical systems, (3d XY modelp, 3d Heisenberg 
modeiy ... )the prototype being the 2d Ising model0,||. Gauge theories are additional 
examples (compact U{1) gauge theoryP,pT|, N = 2 SUSY QGD|n|]), and so are M-string 
theories!^. We will discuss how it works for QGD. 
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The question is: what are the dual excitations of QCD, and what their quantum num¬ 
bers, which determine the symmetry of the conhning vacuum. 

Two main ideas exist: 


a) Monopoles. Their condensation produces dual superconductivity of the vacuum and 
conhnement of electric charges via dual Meissner effect ||T3[| . 

b) Vortices^]. A vortice is a defect associated to a closed line C. If B{C) is the operator 
which creates a vortex on C and W{C') the operator which creates a Wilson loop 
on the contour C", they obey the algebra 


B{C)W{C') = W{C')B{C) exp 


.ncc’r, 


( 6 ) 


where ncc is fhe winding number of the curves C and C. 


/,From the above algebra it follows|^ that whenever (W{C')) obeys the area law 
{B{C)) obeys the perimeter law, and viceversa if {B{C)) obeys the area law, {W{C')) 
obeys the perimeter law. Some people speak of condensation of vortices: however 
the fact that {B{C)) 7 ^ 0 has no special meaning for the symmetry, in the same 
way as the fact that (W{C')) 7 ^ 0. What matters is the area or perimeter law. This 
can be synthetized[|^] in terms of the Polyakov line (L) and the dual Polyakov line 
(t’Hooft line) {L'). Area law for (W) means (L) = 0, perimeter law (L) 7 ^ 0. In the 
same way for {B) area law means {L') = 0, perimeter law {L') 7 ^ 0. 


So for T <Tc one expects (L) = 0 , {V) 7 ^ 0 ; for T > Tc {L') = 0 , (L) 7 ^ 0 . {L') is 
a disorder parameter, (L) an order parameter for conhnement. 


In the presence of quarks however B{C) is not only dependent on the line C, but 
can only be dehned on surfaces having C as a contour, and depends on the choice 
of the surface |1l5|. 


3. Defining monopoles: the abelian projection. 

A conserved magnetic current jff{x) can be associated to any operator (p{x) in the 
adjoint representation. We shall consider for simplicity the case SU{2): all what follows 
can be easily extended to generic SU{N). 

Let us start by dehning (p{x) = (p(x)/|(p(x)|, a direction in colour space, (p is dehned 
everywhere except at sites where 0{x) = 0. According to ref.( 0 ) we can dehne a singlet, 
gauge invariant held strength 

{Df,ip A D^ip)) = ip{df,A^ - - ^(p{di,(p A d^p) (7) 

The two terms are separately gauge invariant, but the coefficients are chosen in such a 
way that bilinear terms in and A^d^p cancel. Going to a gauge p = const (eg. 

p{x) = ( 0 , 0 , 1 )) {abelian projection) the last term in eq.(|^) becomes zero and reduces 
to an abelian held|0 = d^j^A^ — d^A^ The abelian projection is non singular except 
at zeros of p{x) where, after the projection, U{1) Dirac monopoles appear. 
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Defining the dual tensor F*^ = and Jjf = d'^F*^, because of the antisymme¬ 

try of F*^, is conserved 

= 0 ( 8 ) 

In non compact formulations = 0 (Bianchi identities); in compact formulations like 
lattices it can be nonzero |]T7[1. 

The magnetic U{1) symmetry implied by eq.(^ can be 


a) Wigner. The magnetic charge Q is dehned as an hermitian operator and the Hilbert 
space splits into superselected sectors of dehnite magnetic charge. In this case the 
vacuum expectation value (fi) of any operator /i carrying magnetic charge vanishes 
(h) = 0 


b) Higgs broken. Then at least one magnetically charged operator /i exists, such that 
(fi) 7 ^ 0. The free energy density (effective lagrangean) can then be written by 
symmetry and dimensional arguments in terms of {fi) and the dual vector 
potential, as a Higgs Lagrangeanand dual superconductivity follows. 


In conclusion, in non abelian gauge theories, a conserved magnetic current exists, associ¬ 
ated to any choice of the held (p. p is arbitrary, except for the fact that it belongs to the 
adjoint representation. For SU{N) the construction provides N — 1 conserved currents. 

The idea|W,^] is then that the U{1) symmetry is Higgs broken in the conhned phase 
(monopoles condense in the vacuum, electric charges are conhned by dual Meissner ehect), 
becomes Wigner at the deconhning transition and stays such in the deconhned phase. 

The dependence on the choice of p is not a priori known. 

This pattern is detectable by constructing a disorder parameter {fi), the v.e.v. of a 
magnetically charged operator [j^^. Then {fi) 7 ^ 0, T < Tq 



{fi) = 0 T > Tc 


4. Construction of the disorder parameter {fi). 

Expressing p in terms of the helds p of the direct description is an explicit realization 
of the duality transformation. The general idea is to shift the helds by the classical, 
topological conhgurationpll which is created. In the Scrodinger representation | 93 (x)), 
this is done in terms of the conjugate momentum operators n(x), dehned by the canonical 
commutation relations 


[p{xQ,x),Tl{xo,y)] =i5'^{x-y) 

The formula for p{x) is then 

s) = exp (i j(fyn{xo,y)f(y- x) 


(B) 

( 10 ) 


giving 


p(x)jp(xo,x)) = \p{xq,x) + p{x)) 


( 11 ) 
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Eq. (^) is the analog of the familiar translation for a 1 dimension system 
= \x + a) 


( 12 ) 


The construction of fi for monopoles in compact U{1) gauge theories, in SU{N) gauge 
theories p]| , P^ , in full QCD[^ and for vortices in QCD|^ has been developed in a series 
of papers. Similar constructions have been done for spin system exhibiting duality 
Most of these systems have well known phase structure, and the construction was made 
to test the method sketched above, a test which proved successful! in all cases. We will 
list here the problems encountered and solved, and then present the results obtained and 
discuss their implications. 

1. Compactness In compact theories the helds are dehned in a hnite range, and cannot 
be shifted by arbitrary amount, as in eq. (pT]) . 

We shall discuss compact U{1), with Wilson’s action 

S'=/3 ^ (1 - cos6'^i.(n)) Z = exp(-S') 

The creation operator of a monopole at x and time xq, can be written in this case 


fi{x) = exp 


{cos(6'o*(n,xo) - bi{x - n)) - cos(6'o*(n, xq)} 


(13) 


For small angles 6oi,bi -C 1, /i(x) ~ exp[/36oi{n,xo)bi{x - n)], 6^^ = is 

the held strength tensor, = Ei the electric held, i.e. /i(x) reduces to the usual non 
compact formula (pT|) for small helds. When computing (/x) = /['D6*]/iexp(—S)/Z it is 
easily realized that 


{,■) = ZjZ 


(14) 


where Z = exp(—S) and S is obtained from S replacing cos(6'oi(fi, xq)) at time xq in the 
action by cos(6'oi(fi, xq) — bi{x — n)). 

It is a theorem |]I0| that fi creates on the lattice a monopole at t = xq which then 
propagates at f > Xq. The form (|^ for the disorder parameter is common to all systems 
with duality. 

2. Fluctuations {fi) = ZjZ is not easy to compute numerically like any partition 
function. Indeed {fi) is (eq ([T3|) ) the vev of the exponential of an integral extending to the 
3d space, which is roughly proportional to iV|, and huctuates as producing for {fi) 

huctuations C>(exp(iV|'''^). The way outp^ is to compute the quantity 




(16) 


which is the diherence of the expectation values of two actions. All the relevant infor¬ 
mation contained in {fi) can be extracted from p. The thermodynamical limit is reached 
by hnite size scaling analysis. A typical shape of {p) is shown in £g.(l), which describes 
the deconhning transition in compact f/(l)[|lOl. {p) is 0{i) for jS < I3c and drops by 50 
orders of magnitude at jSp in the inhnite volume limit the fall becomes sharper and {p) 
is 0 for T > Tc- The corresponding form of is shown in £g.(2). 
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Figure 2. Fig.2 p for compact U{1) 


Figure 1. Fig.l (p) for compact U{1). 


In terms of p (p) = exp Jq p{x)dx . For P < P goes to a finite limit as V oo. 
For P > Pc, p ^ —oo, as 1/ —>• cx). Indeed p = —kNs + k', with k > 0 and Ns the 
spatial extension of the lattice: this corresponds to hmy^oo(/w) = 0, P > Pc- 
Around Pc 



(16) 


6 being the critical index of the order parameter. 

3. Locality In principle the sum at the exponent in eq.(p!3D could produce infrared 
problems, e.g. destroying cluster property 


~ Cexp(-/l|a:|) + |{p)p 


( 17 ) 


which is needed to dehne properly (p). 

4. Gauge invariance The vev of a (magnetically) charged operator can have problems 
with gauge invariance [Q . 


All these problems have been solved and checked on the systems where the result is 
known, see eg. 0. 

We will briefly comment on point 3 and 4. 


5. Gauge invariance 

In a gauge invariant formulation of a U{1) gauge theory, the ground state is gauge 
invariant (17a|0) = |0)). If (p{x) is a charged local operator Ul^ip{x)UA = ip{x) exp{iA{x), 
and hence 

{0|43(x)|0) = {0 |£/P(x)C/a| 0) = e“<')(0|^(n|0) VA 


(18) 
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Therefore (0|(^(a;)|0) = 0. 

The lattice version of this statement is known as Elitzur’s theorem||24||. 

The usual way out is to fix the gauge (e.g. the unitary gauge) and look for spontaneous 
breaking of U{1) symmetry ||26||. A lattice formulation is per se gauge invariant. How can 
then a Higgs phenomenon be described? 

The solution consists in dehning a gauge invariant charged operator (p{x), and in taking 
(0|(^(x)|0) as an order parameter. The dehnition of cp is the following 


(p{x) = ip{x) exp[z(h, A)] 


(19) 


where 


{h,A)=JdSA^,{y)h^{y-x) and d^hf,{z) = 5^{z) 


Under a gauge transformation A(a;), with A (a;) —0 as |a;| —> cx) > A^ + A and, 
by partial integration (h, A) —> (h, A) — A(a;), while p{x) —> The phase A(a;) 

cancels between the two factors of eq. ([I^ and (p{x) —> (p{x)\ (p{x) is gauge invariant. 


Under a global transformation instead A^ 

~ i A ~ 


A^, (h, A) — {h,A), (p —>• and 




( 20 ) 


The exponential factor in eq. m makes p non local. 

A popular choice for corresponds to a Schwinger parallel transport. E.g. if the string 
is put along the time axis h^(z) = 6q9{zo)6^{z) and 


Ps{x) = p{x)exp (i f A^{x,T)d 

\ J OO 

Another possibility is Dirac choice [P7 


( 21 ) 


A^iz) = {0,Hzo)C(z)) 
for which 


C{z) = 


dvr |T|" 


Pd{x) = p{x) exp W d^yA{y,xo) 


1 y — X ' 
dvr 1^— 


( 22 ) 


It can be shown that the correlator ((p5(x)(p5(0)) has infrared problems, in that it van¬ 
ishes exponentially at large distances {(ps{x)(psiS^)) ~*’|a;|^oo Aexp(—p|a:|) even in the Higgs 
phase, due to the presence of the string. It is then not suited to dehne the order parameter. 
The Dirac choice instead is local enough to allow cluster property, and 


(^d(2:)^d(0)) + 2 lexp(-p|x|) 


(23) 

defines the order 


with {(p) 7 ^ 0 in the Higgs phase, {(p) = 0 in the Coulomb phase. Eq. 
parameter. 

The operator p which creates a monopole ||10|| is gauge invariant,magnetically charged, 
and Dirac like [|], 3^ 


No infrared problem nor problem with gauge invariance exists. 
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Figure 3. Fig.3 p for monopole condensation 

in SU{3). Figure 4. Fig.4 Finite size scaling of p. 


6. Results for confinement. 

We can summarize the results as follows. For quenched QCD (p) can be defined by a 
recipe analogous to that of compact in any abelian projection and p = -^ 

can be determined. The typical behaviour of p is shown in hg.(3) and is qualitatively 
similar to that in hg.(2) for 17(1). 

In order to perform the inhnite volume limit p is measured for different spatial sizes of 
the lattice. 

For f3 < f3c, p goes to a hnite limit as Ns —> oo. 

For P > Pc, P ^ —kNs — k' {k > 0), i.e. Umjsfg^ooip) = 6xp(/^ p{x)dx) = 0 
At /3 ~ Pc, by dimensional arguments, 
r Nc CL Nt 

where ^ ^ {P — Pc)~’^ is the correlation length, which goes large with the index u as 
P —> Pc- Near Pc a/^ ~ 0, Nt/Ns ~ 0 and {p) ~ T^^{Ns^^r) or 


P 


= f + HN'/t) 


r 

This scaling law is obeyed, fig. (5), with the known values of Pc and z/ and allows to 
determine S. 

The result for 517(3) is 5 = 0.50(2) independent of the abelian projection. 

The computation has been performed in a number of abelian projections, and also 
by choosing ||23||, conhguration by conhguration, as direction of abelian projected $ the 
conventional 3 axis used to prepare the conhguration. This gives an average on a number 
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of abelian projections eqnal to the nnniber of configurations used, and again gives the 
same behaviour for (/i). 

A similar operator p can be defined which creates a t’Hooft line. Fig. (6) shows a 
comparison between the corresponding p and that of monopoles, showing that they almost 
overlap. 


7. Full QCD 


In full QCD (p) can be defined in the same way as in quenched, by modifying only 
the gauge part of the action in the definition of S, eq.(|^). Simulations give a behaviour 
similar to the quenched case, and a negative peak at a /3 which coincides with that of the 
chiral transitions^ . 

A full hnite size scaling analysis is on the way, as well as an analysis of p at S> /9c to 
establish that (/i) is a correct disorder parameter. 

In principle dual j ^ :.li_ .li quarks, so 

that (p) would be < ). 



Fig.5 p for monopole condensation (full cir¬ 
cles), and chiral condensate (open circles) for 
2 flavour QCD. 


8. Conclusions and outlook. 


1) deconfinement is an order disorder transition, and the mechanism is dual supercon¬ 
ductivity. Monopoles condense in the confined phase in all abelian projections, cer¬ 
tainly in quenched QCD, and also, by preliminary evidence, in full QCD, Nq —> C )0 
limit is respected. 

2) The dual configurations of QCD are not known. Whatever the dual configurations 
are, they carry magnetic charge in all abelian projections. 


3) The dual Polyakov line is a good disorder parameter, and has the same behaviour 
as the monopole (p) |T^ . 


The analysis of full QCD must be completed.lt will require much computer time and 
work. 
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A breaktrough is needed to identify the dual excitations. A possible direction could be 
that indicated in ref.[^, in which a disorder parmeter is defined which coincides with 
ours 0{a^) in all abelian projections. 

I am grateful to L. Del Debbio, M. D’Elia, B. Lucini, G. Paffuti for discussions and 
help in preparing this talk, which is based on our common work. 
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